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The Anderson model of a twofold spin degenerate impurity level in the limit of infinite Coulomb 
repulsion, U ^ oo, coupled to one and two degenerate conduction bands or channels, is considered 
in pseudo-particle representation. We extend the Conserving T-Matrix Approximation (CTMA), 
a general diagrammatic approximation scheme based on a fully renormalized computation of two- 
particle vertex functions in the spin and in the charge channel, to the calculation of thermodynamic 
and spectral properties. In the single-channel case, the CTMA yields in the Kondo regime a tem- 
perature independent Pauli spin susceptibility for temperatures below the Kondo temperature Tk 
and down to the lowest temperatures considered, reproducing the exact spin screening in the Fermi 
liquid state. The impurity spectral density appears to remain non-singular down to the lowest 
temperatures, in agreement with Fermi liquid behavior. However, the unitarity sum rule, which 
is crucial for an impurity solver like the CTMA to be applicable within Dynamical Mean Field 
Theories for strongly correlated lattice models, is overestimated at the lowest temperatures. We ar- 
gue that this shortcoming may be due to numerical imprecision and discuss an appropriate scheme 
for its correction. In the two-channel case, the spectral density calculated within CTMA exhibits 
qualitatively the correct non-Fermi liquid behavior at low temperatures, i.e. a powerlaw singularity. 

PACS numbers: 71.10.-w, 71.27.4-a, 71.28.-|-d, 71. 55.-1 
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I. INTRODUCTION 

Over the past two decades the problem of correlated 
electrons on a lattice has emerged as a central theme of 
condensed matter theory. With the exception of one- 
dimensional systems, there are no systematic analytical 
methods available for solving models like the Hubbard 
model. A powerful approximation scheme is the Dynam- 
ical Mean Field Theory (DMFT), in which the lattice 
problem is mapped onto an effective single impurity An- 
derson model (SIAM), with self-consistently determined 
properties of the conduction electronsi2(. It is a non- 
trivial task to solve these SIAMs. The properties of 
real Kondo or mixed valence impurities in metals are 
of interest in their own right, with recent emphasis on 
non-Fermi liquid behavior in multi-channel modela^i^i^. 
Various methods have been successfully applied to solve 
these models in certain parameter regimes. The Bethe 
ansatz (BA) method allows to calculate the thermody- 
namic properties of models with a flat conduction elec- 
tron density of states^. Bosonization methods have 
been used to obtain, for example, the finite-size spec- 
trum of one- and two-channel impurities^. Conformal 
field theory is a powerful tool to analyze the low energy 
excitations of multi-channel modelsiSiiLiS. The method 
of continuous unitary transformations"'^' has been suc- 
cessfully applied to the Kondo model in the vicinity of 
the Toulouse poin^. These analytical methods are com- 
plemented by numerical methods like Quantum Monte 



Carlo (QMC) simulations (for not too low temperatures 
and moderate U)i^*iSiii and Wilson's Numerical Renor- 
malization Group (NRG) which has been very success- 
ful for not too large degeneracies in the spin or charge 
channeli2*i2i2a. 

The difficulty with quantum systems of the Anderson 
impurity type is the strong on-site Hubbard repulsion, 
which effectively constrains the quantum dynamics to 
a Hilbert space with fixed impurity occupation num- 
ber and makes these problems inaccessible by straight- 
forward perturbation theory. It is, in particular, difficult 
to describe the weak coupling (fluctuating local moment) 
behavior at high energies and the strong coupling fixed 
point behavior, realized below a strong coupling energy 
scale, typically the Kondo temperature Tk, by a single 
technique. In view of possible applications as an "impu- 
rity solver" within DMFT methods or to quantum im- 
purity and quantum dot systems with a complex local 
spectrum, an accurate method which does not rely on 
integrability conditions or on the simplicity of the local 
or conduction electron spectrum is highly desirable. 
For that purpose we had proposed earlier a general 
diagrammatic approximation scheme'^-. The starting 
point is a pseudoparticle representation of the impu- 
rity level, where the constrained dynamics are built into 
the very definition of the quantum fields22,, and approx- 
imations conserving its internal symmetry are defined 
by means of Luttinger-Ward functionals. The conserv- 
ing approximation which incorporates the dominant, lo- 
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cal spin and charge fluctuations on the level of a fully 
renormalized calculation of the total two-particle ver- 
tex (or T-matrix) has been termed the Conserving T- 
Matrix Approximation (CTMA). In contrast to previ- 
ous approximations like the Non- Crossing Approxima- 
tion (NCA)SiSiSi and its extensiong^"?:, the CTMA de- 
scribes the weak as well as the strong coupling behavior 
of the single-channel SIAM correctly on the level of the 
pseudoparticle propagators2ii2i. However, physically ob- 
servable spectral properties had not yet been calculated 
because of their computational complexity. 
In this article we present CTMA results for the thermo- 
dynamic and spectral properties of the SU(A^)xSU(M) 
Anderson impurity model, N being the local spin de- 
generacy and M the number of identical, conserved con- 
duction channels. We will focus on the single-channel 
Fermi liquid case {N = 2, M — 1), although results 
for the two-channel non-Fermi liquid sector of the model 
(TV = 2, M = 2) will also be shown. The spin susceptibil- 
ity as well as the frequency dependence of the impurity 
electron selfenergy indicate that the spin-screened Fermi 
liquid ground state of the TV = 2, M = 1 SIAM is in- 
deed captured by CTMA. However, the unitarity sum 
rule of the spectral density, which is vital for DMFT 
applications, is overestimated. A detailed inspection of 
the impurity electron selfenergy shows that this failure 
seems to originate from an imprecise treatment of high- 
energy processes, either due to numerical inaccuracy or 
due to CTMA neglecting non-singular potential scatter- 
ing terms, and that such imprecision influences the low- 
energy behavior via the Kramers-Kronig relation. Based 
on this analysis we propose below a phenomenological 
correction scheme which imposes the causality of the im- 
purity selfenergy, and which may thus make the CTMA 
applicable as an impurity solver for DMFT calculations. 
This approach amounts to adding an appropriate poten- 
tial scattering term to the real part of the impurity self- 
energy, taken to be a temperature independent constant. 
We will term this scheme the "effective potential scatter- 
ing method" . 

The paper is organized as follows. We describe the con- 
serving pseudoparticle technique in section ^1 includ- 
ing several physical and technical justifications of the 
CTMA. The detailed CTMA self-consistent equations are 
given in the Appendix. The CTMA results for the tem- 
perature dependent, static spin susceptibility and for the 
impurity spectrum are presented in sections IIIII and IIVI 
respectively. Our effective potential scattering correction 
scheme is discussed in detail at the end of section Hvl In 
section O we compare the CTMA and NCA result for 
the spectral density of the two-channel SIAM, where the 
ground state is not a Fermi liquid. We conclude with a 
discussion of the results in section IVII 



II. MODEL AND CONSERVING T-MATRIX 
APPROXIMATION 

We consider the SV{N)xSV{M) SIAM in the hmit of 
infinite Coulomb repulsion, implying that the Af-fold de- 
generate impurity level (called d- level here), labeled by 
spin a = —N/2, . . . ,+N/2, is at most singly occupied. 
The empty impurity state is M-fold degenerate, labeled 
by fi — 1, . . . , Af, and is coupled to a corresponding de- 
generate degree of freedom in the conduction band, e.g. 
z-component of angular momentum. In the pseudopar- 
ticle representation, the singly occupied (empty) level is 
created by fermionic (bosonic) operators (&J,), which 
satisfy the constraint Q ^ Y.a fUa + T^^b^b^ = 1. 
The physical (or d) electron creation operator on the im- 
purity site is 4 = Y^aflh- The SU(iV)xSU(M) An- 
derson impurity Hamiftonian is then defined by 

kcr/x CF 

+ Vj2i^l^by^ + h.c.) + XQ. (1) 

Here, Cq^^ = X^k '^ko-^ creates a conduction electron at 
the impurity site R — and e^.o- = + o'gfJ'sB is the 
impurity level in a magnetic field with fj,B and g the 
Bohr magneton and the Lande factor, respectively. The 
operator transforms according to the conjugate rep- 
resentation of SU(M). We denote the density of states 
(DOS) of the conduction electrons at the Fermi energy ep 
by iV(0) and assume it to be structureless. (All numerical 
results were obtained for a Gaussian DOS. Note, how- 
ever, that our method works for arbitrarily structured 
DOS.) Instead of the hybridization V, we will frequently 
use r — 7r7V(0)y^ as a parameter of the model. In the 
Kondo regime, F <C where the low energy excitations 
resemble those of the Kondo model, the above Hamilto- 
nian has a dynamically generated strong coupling scale, 
the Kondo temperature Tr-, where perturbation theory 
breaks down, 

^^' = ''l^J ^"pI-atJ' (2) 

where 2D is the bandwidth of the conduction electron 
DOS. The Hamiltonian ^ possesses a U(l) gauge sym- 
metry with respect to simultaneous transformations of 
/ct and bp related to the conserved charge Q, A being the 
local gauge field. The exact projection of the dynam- 
ics onto the Hilbert subspace Q = 1 is accomplished by 
taking a gauge with a time independent A, re-defining 
the zero of the energy scale as uj uj + X, and letting 
A — > oo in all expressions; see Ref. 28] and the appendix 
of Ref. for details of the projection technique. The 
charge conservation in conjunction with the constraint 
Q = 1 implies an orthogonality catastrophe between the 
Q — 1 initial and the Q — final states, and leads to 
infrared threshold power-law behavior of the pseudopar- 
ticle Green's functions Gf^b{^) <x w""-'-''. In the Fermi 
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liquid case, M < N — 1, the exponents a/^f, are closely 
related to the average impurity occupation number n^i 
via the Friedel sum rul o^^i??'?^ , 



(3) 



Approximation schemes for calculating Gf^ which vio- 
late the gauge symmetry would, hence, violate the or- 
thogonality of initial and final states, and should be 
expected to give incorrect results for physical quan- 
tities, even though certain aspects of the Fermi liq- 
uid fixed point can be described by symmetry-breaking 
approximation a?^i??i?'*i?^'??i?7i?? . Therefore, we take 
great care to preserve the gauge symmetry. It can be 
reconciled with the time independent choice of the gauge 
field A necessary for the Q = 1 projection by employing 
a conserving approximation, derived from a Luttinger- 
Ward generating functional The local selfenergies 

Set, a = f,b, c, defined by 



G-J{oj) = 



(4) 

(5) 

(6) 



where Gj^ 



— X]fe(^ ^ ^fc)^^' ^'"S generated by 
functional differentiation of $ with respect to the self- 
consistently renormalized Green's functions Y^aiw) = 
6^/6Ga{i-^)- As a result, for any given approximation 
to $ a set of self-consistent nonlinear integral equations 
for and Ef, is obtained, which, in general, cannot be 
solved analytically but is amenable to numerical solu- 
tion. The central task is then to find the correct gener- 
ating functional which captures the essential physics of 
the problem at hand. 

Non-Crossing Approximation (NCA). The NCA is often 
used for its computational simplicity to obtain a rough 
description, and it even captures the universal behavior 
inherent to Kondo-type problems^. However, the NCA 
recovers the correct Kondo scale Tk only because of a 
fortunate compensation of the neglect of spin-flipping 
logarithmic terms and an incorrect logarithmic resum- 
mation of potential scattering termsi2i_''^'^ . Below Tk it 
develops spurious infrared singularities in physical quan- 
tities. The NCA breaks down in a magnetic field B even 
in the weak coupling regime (T > Tk and/or B > Tk), 
producing, in addition to the two Zeeman-split Kondo 
peaks a third, field independent resonance in the im- 
purity spectral density at the Fermi energy. It seems 
that this spurious behavior originates from the incorrect 
treatment of the potential scattering in NCA mentioned 
above. On the level of auxiliary particles, the NCA does 
not give the correct FL threshold exponents Eq. ^ , but 
instead af^^ = M/{N + M), a^^^ = N/{N + M). 
It can be shown by power counting arguments, using 



,NCA 



= 1, that any self-consistent calculation 
involving only a finite number of skeleton selfenergy dia- 
grams just reproduces the incorrect NCA exponents. 
Conserving T-Matrix Approximation ( CTMA ). Hence, 



selfenergies and two-particle vertex functions must be 
comprised of an infinite class of skeleton diagrams in or- 
der to describe the Fermi liquid fixed point. Since the 
latter is a consequence of the singlet formation between 
conduction electron and impurity spin, it is natural to 
assume that higher than two-particle correlations need 
not be considered in the single-channel case. The total 
vertex functions of conduction electrons (c) and local de- 
grees of freedom (pseudofermions /, slave bosons b) are 
then two-particle T-matrices which are obtained from 
an infinite summation of irreducible parts via Bethe- 
Salpeter equations. We use the smallness of the pa- 
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FIG. 1; Diagrammatic representation of the Bethe-Salpeter 
equations for the vertex functions t''^^' and t''^''' defining the 
CTMA. Dashed, wiggly, and solid lines represent here and in 
the following the renormalized pseudofermion, slave boson, 
and the local conduction electron propagators, respectively. 
The external lines are for clarity only and are not part of the 
vertices. 

rameter VN(0)^ 1 to select the leading diagrams of the 
irreducible parts. This results in the ladder approxima- 
tion for the total two-particle vertices shown in Fig. ^ 
The Luttinger-Ward functional that generates by second 
functional differentiation the vertex functions of Fig. ^ 
is constructed by connecting the entry and exit points 
by Green's function lines and is shown diagrammatically 
in Fig. 121 The diagram containing two (renormalized) 
boson lines is not a skeleton, is already contained in the 
first (NCA) diagram via self-consistency, and, hence, is 
omitted. The conserving approximation obtained in this 
way has been called Conserving T-Matrix Approxima- 
tion (CTMA). The self-consistent equations for the ver- 
tex functions and selfenergies to be solved are given ex- 
plicitly in Appendix fXl Note that the / — c and the b — c 
vertices in Figs, ^andl^ldescribe spin and charge fluctua- 
tions, respectively. Therefore the CTMA should provide 
a good approximation not only in the Kondo, but also in 
the mixed valence and empty impurity regimes. 
On a more formal level, the CTMA can be justified both 
near the weak and near the strong coupling fixed points. 
Expanding, in the weak coupling regime, the CTMA 
in terms of hare, projected Green's functions {B — 0), 
G'O^(w) = l/(w±iO), Gl{uj) = I/(w + ed±iO), it is seen 
that the CTMA c-f vertex is exact up to leading loga- 
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coupling fixed point of the single-channel SIAM as well. 
This will be shown by numerical evaluation in the follow- 
ing sections. 
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FIG. 3: Diagrams containing the leading logarithmic contri- 
butions to the c-f vertex function. In contrast to Figs. Q |5] 
the lines represent the bare Green's functions here, a) Con- 
tribution from the NCA generating functional, b) additional 
contribution from the CTMA. Spurious logarithmic potential 
scattering terms cancel only when both terms are added. The 
c-b vertex has no logarithmic terms and, hence, does not flow 
under perturbative RG. 
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FIG. 2: Generating functional of the CTMA. The first dia- 
gram corresponds to the NCA generating functional. Each 
diagram contains exactly one closed auxiliary particle loop; 
diagrams with more than one such loop vanish after projec- 
tion onto Q = 1. The class of CTMA diagrams beyond NCA 
is uniquely defined by the condition that each conduction 
electron line spans exactly two bare three-point hybridiza- 
tion vertices, as shown. The diagrams of the first and of 
the second column generate, by second functional derivative, 
S'^^/SGcSGf and S^^/SGcSGt, the spin and the charge fluc- 
tuation T-matrices, T*"^^', T<'=''\ respectively (Fig.[Tll. 



rithmic order, as seen in Fig. |31 Therefore, CTMA does 
incorporate the correct renormalization group (RG) flow 
in the weak coupling region^^. In particular, logarith- 
mic potential scattering terms, present in each one of the 
diagrams of Fig. 13 cancel correctly when the two dia- 
grams are added. Therefore, we expect that the CTMA 
correctly describes the Zeeman splitting of the Kondo res- 
onance even in a large magnetic field, in contrast to the 
NCA, where only the first diagram of Fig. |31is included, 
and its logarithmic potential scattering part leads to a 
spurious third peak at the the Fermi energy oj = (see 
above). 

By virtue of the self-consistent inclusion of selfenergy di- 
agrams in the propagators of Fig.^the CTMA vertex not 
only includes ladder but also parquet diagrams. More- 
over, at any given order of self- consistent perturbation 
theory, the CTMA includes all diagrams with the max- 
imum number of local spin and charge fluctuation pro- 
cesses in the sense of principal diagrams. Since these are 
expected to be responsible for the formation of the spin- 
singlet in the Kondo or mixed valence regime, CTMA 
may be expected to capture the physics of the strong 



III. SPIN SUSCEPTIBILITY 

In this section we report the CTMA result for the static 
impurity susceptibility 



(7) 



5=0 



of the A = 2, A/ = 1 SIAM in the Kondo regime and 
compare with NCA and Bethe ansatz results. Through- 
out the paper, the evaluations in this regime will be done 
for the typical parameter set ed/D = —0.81, T/D — 0.2. 
In Eq. yielding a low-temperature impurity occupa- 
tion number of n^a ^ 0.47 and Tk = 4.16 x lO'^D. 
Other parameter sets in the Kondo regime give similar 
results. M = gfJ-s'^a^ o"n^ is the impurity magnetiza- 
tion. 
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(8) 

is the impurity occupation number with spin a projected 
onto Q = l.^&^ (/? = l/fcsT), and the magnetic field B 
couples only to the impurity spin (Eq. Q). The expres- 
sion {Tj) is equivalent to the w = limit of the causal 
dynamical linear response susceptibility 

X^{T,UJ = 0) = -iJ dteit){[M{t),MiO)]) . (9) 

This is readily shown from Eqs. |7J, employing 
d/dB = [dGf^/dB) S/SGf^ and 
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dB 
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^ SGf^. dB 
which follows from the definition of G/o-, 
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is the irreducible four-point pseudofermion vertex. Any 
conserving approximation by construction fulfills the 
equivalence of Eqs. Q and since it respects the Ward 
identity . We choose to use Eq. Q , because it is com- 
putationally less demanding than the correlation function 
Eq. ®. 

Note that additional terms Xb{T) and Xib{T) arise, if B 
couples also to the conduction electron spin. Xb{T) is 
the constant Pauli susceptibility of the conduction band 
and XibiT) a mixing term correlating the impurity and 
the conduction electron magnetization. Since the latter 
is, for a flat DOS, due to the electronic polarization at 
the bottom of the band, Xib{T) is usually negligible for 

Xi(r) is of principal interest as an indicator of whether 
CTMA captures the spin singlet Fermi liquid ground 
state of the single-channel SIAM. The result is shown 
in Fig. 0] While at exponentially high temperature, 
ln(T/TA-) » 1, x^{T) = (l/4)5VI/r, typical for a free, 
fluctuating local moment, Xi{T) shows T- independent 
Pauli behavior for T < 0.1 Tk and down to the low- 
est T considered, characteristic for the completely spin- 
screened Fermi liquid state. By contrast, the low-T be- 
havior of the NC A is a power law, xf^^ (T) - xf^^ (0) oc 
At T = 0, Xi{T) acquires the value 

X40) = , (12) 

which defines the universal low-temperature scale, T^, 
of the Kondo or Anderson model, related to Tk by the 
Wilson number W = Tr/Tl- 

The comparison of the CTMA result with exact 
methods like NRCii^ or BAiS*^ can be made quan- 
titative. The dimensionless quantity XiC^/^-R") = 
Xi(T)/{g'^^\/TK) is known to be a universal function 
ofT/Tx (i.e. independent of the microscopic parameters 
of the Hamiltonian) , with Xj(0) = W. Self-consistent 
approaches like NCA or CTMA reproduce this universal- 
ity, since they include a resummation of the logarithmic 
terms of perturbation theorji^. In comparing the CTMA 
and the BA results one must, however, observe, that the 
breakdown scale of perturbation theory, Tk, depends on 
its precise definition. Therefore, care must be taken that 
the same definition of Tk is used for both, the CTMA 
and the exact method. In Wilson's original work on the 
Kondo modeli^, a Kondo temperature was defined 
such that in the high temperature expansion of Xi (T) all 
terms of 0(ln(r/T^)^^) cancel each other. Rasul and 
jjg.^gpj-^46j47j48 ygg(j ^}ig same criterion for the SIAM and 
found for the Kondo temperature, 

where C = 0.5772157 is Euler's constant. With this def- 
inition, the universal Wilson number was found to be 
W=0.4128. Using the same definition, we find within 



CTMA (Fig. mi, 

^(CTMA) ^ Tl,xm/t^l = 0.462 . (14) 

In the BA method, Tk is defined in a somewhat differ- 
ent way due to a different cutoff scheme, resulting in a 
Wilson number W^^ = exp{C + l/4)/7r3/2 w 0.4107^''*4 
(for N ^ 2). Therefore, in the BA curve for x^{T/Tk) 
we rescale Tk such that Xi(0) obtains Wilson's T = 
value W — 0.4128 (Fig. As seen from the figure, 
the CTMA result for the static susceptibility is in strik- 
ingly good quantitative agreement with the BA result 
not only for T > T^, but also in the strong coupling re- 
gion, r < 0.1 T^. This shows that the CTMA describes 
the low-energy excitations around the Fermi liquid fixed 
point even quantitatively correctly at least in a thermo- 
dynamic quantity like the magnetic susceptibility. From 
the general properties of conserving approximations, one 
may expect the same to be true for dynamical quantities 
as well. This will be investigated in the next section. 
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FIG. 4: The static magnetic susceptibility as a function of 
temperature. The temperature and the susceptibility scale 
are plotted in units of the conventional Kondo temperature 
Eq. (Uni, that leads to the value W = 0.4128 for the 
Wilson number, see text. The CTMA curve agrees well with 
the exact BA result in both, the high temperature and in the 
strong coupling regimes, with only a small spurious interme- 
diate maximum in the crossover region. 



IV. DYNAMICAL PROPERTIES 

The Pauli behavior of the impurity susceptibility cal- 
culated in the previous section shows that the spin struc- 
ture of the low-T excitations captured by the CTMA free 
energy functional is such that it describes the complete 
screening of the impurity spin by the conduction elec- 
trons correctly even with good quantitative precision. It 
is the resulting absence of spin fiip scattering at energies 
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below Tk that is responsible for the Fermi liquid behavior 
near the strong coupling fixed point of the single-channel 
Kondo or Anderson model. Therefore, one may conjec- 
ture, that the CTMA may capture the Fermi liquid na- 
ture of the low-T excitations as well, with well defined 
quasiparticles which should become visible in dynamical 
quantities like the impurity spectral density and the self- 
energy. This expectation is further supported by the fact 
that the CTMA does indeed reproduce the correct Fermi 
liquid threshold exponents on the level of the auxiliary 
particle dynamicsSi. 

The quantities of prime interest, e.g., for DMFT are 
the impurity electron Green's function. 



Lu - ed±iT - ± lO) 



(15) 



and the interaction part of the selfenergy, Y,da{'^)i in the 
spin-screened case, N = 2, M — 1. The Fermi liq- 
uid theory implies certain exact low-energy properties, 
which the CTMA results must be compared to, namely 
the unitarity limit for the impurity spectral density at the 
Fermi level A(i„{Q), the half width T and position l^. of 
the Kondo resonance, as well as the low-energy behavior 

of ln\Y.dcr{l-^), 
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where i? = 2 is the Wilson ratio. The proof of these 
relations is compiled in Appendix IbI 



A. CTMA solution and Fermi liquid behavior 

For the computation of Gda{^) within the CTMA, ob- 
serve that it is related to the single-particle conduction 
electron T-matrix, tccr('^) = V^Gdaii-^), where, after pro- 
jection onto Q ~ 1 only diagrams with a single pseu- 
doparticle loop, i.e. irreducible diagrams, remain. In the 
conserving scheme it is, therefore, constructed as 



Gdai(^) 



1 



lim 



1 



(5$ 



y2 A-*oo Q(A) SGcaiuj) 



(21) 



The corresponding CTMA diagrams are shown in Fig. 
|S1 and the details about their evaluation are given in 
Appendix ^ 

It is worth noting that, by definition, the im- 
purity electron propagator is equivalent to the 



/ — &- "particle- hole" correlation function, Gda{t) = 
-i{f{b^t)f^{t)p^{0)b (0)}). This might seem to of- 
fer another possibility of calculating Gda- using the irre- 
ducible f — b vertex Tfb = 5"^^ / 8G jaSGb- However, any 
diagram of the f — b particle-hole correlator constructed 
in this way contains necessarily two pseudoparticle loops, 
and, hence, vanishes by projection. Therefore, in CTMA 
the f — b correlation function is comprised of the (NCA- 
like) f — b bubble diagram only, which is clearly not suf- 
ficient to recover the Fermi liquid strong coupling fixed 
point. We note in passing that the non-trivial contribu- 
tions to Tfb comprising the full Gd are generated from 
free energy diagrams which contain more than one pseu- 
doparticle loop, but are not contained in CTMA. As a 
conclusion, in CTMA Gdai^) must be calculated unam- 
biguously using Eq. (|21|l . 

The CTMA results for the d-electron spectral func- 
tion, Adaii-^) — TaiGdai^ ^ iO)/TT are shown in Figs. El 
anddfor T ^ Tk and T O.OITk, respectively, together 
with its decomposition into the f ~ b bubble contribu- 
tion (1st diagram in Fig. and the vertex corrections 
(2nd and 3rd diagrams in Fig. |SJ| . Notably, even at el- 
evated temperature, T = Tk, the vertex corrections are 
not negligible. For T ^ the f — b bubble develops 
an infrared powerlaw divergence oc |a;|"^"'""''~^. Since 
for our parameter set in the Kondo regime the exponent 
a f + ab — 1 =~ 0.056 is rather small, and the pseudopar- 
ticle exponents, Eq. |21 are obtained asymptotically, this 
singularity starts to develop for T = O.OIT^ merely as a 
discontinuity in the slope of Adaii^)- Most importantly, 
however, the total d— electron spectral function does not 
show any infrared singularity on the scale of T, and the 
width F and position of the Kondo resonance are in 
excellent agreement with the Fermi liquid predictions, 
Eqs. H17|) . H18|) . given the uncertainty in these quantities 
arising from the fact that the Kondo resonance deviates 
from the Lorentzian shape for lu > Tk (see Appendix 
IbI . In contrast, the unitarity limit, Eq. H16|l . is signif- 
icantly violated in Fig. [7| To investigate the origin of 
this failure, we plot in Fig. |S1 the imaginary part of 
the interaction selfenergy I^dai'jJ — *0). It is seen that 
even at the lowest temperature the CTMA result does 
not develop any singularity, in contrast to NCA. How- 
ever, the position of the minimum of ImSdai^ — iO) is 
incorrectly shifted to a negative frequency ujq of 0{Tk), 
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FIG. 5: Diagrams defining the d-electron Green's function 
within CTMA. For details of their evaluation see Appendix 

El 



7 



2 - 



^- 




FIG. 6: CTMA and NCA results for the local spectral den- 
sity Adtriu}) for ea/D = -0.81, F/D = 0.2 at T = Tk = 
6.16 X 10"^ D. The decomposition of the CTMA result into 
/ — 6 bubble and vertex corrections arising from T^-^-f^ (fc ver- 
tex) and t''^''' (be vertex) is also shown. The NCA shows an 
incorrectly large shift of the Kondo peak of 0{Tk) towards 
positive frequencies, which is corrected by CTMA; see the 
text. 



[uj — uj^f behavior for |w — wo| ^ Tk (Fig.lHJ. Its prefac- 
tor, determined from our parameter set from Fig. |3 as 
0.0244 is in excellent agreement with the exact Fermi 
Uquid value, aV = 0.0239 D (Eqs. UHl, The tem- 

perature dependence of the minimum of Iml]dcr(w — iO) 
is analyzed in Fig.^l Again, the CTMA solution shows 
behavior from the lowest T considered up to T ~ Tk, 
ImEdcr(wo — iO) = aT{'nT /TkY , where the prefactor 
aF — 0.013 D is of the same order of the exact value 
ciF, although roughly a factor of 2 too small. 
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FIG. 7: CTMA impurity spectral function and decomposi- 
tion into / — fe bubble and vertex contributions for the same 
parameters as in Fig. |S| but at T = O.OITk. The exact val- 
ues for AdCT(O), f. Id, Eq. CHJ, CJl, CHJ, are displayed for 
comparison. 



where ImEdcr(wo — iO) acquires a spurious negative value. 
Even at the lowest T considered, ImE(jo-(w — VS) shows 



FIG. 8: Correlation part of the physical d— electron selfen- 
ergy Im Etjo- (lo) calculated in CTMA for various temperatures 
and ed/D = —0.81, T/D = 0.2. For the lowest temperature 
shown, r = 0.16 T/f, the NCA resuh (dot-dashed line) is 
shown for comparison, showing the cusp- like infrared singu- 
larity typical for NCA. The bold solid line centered at the 
origin represents the T = behavior expected for our param- 
eter set based on Fermi liquid theory, Eqs. 1191 . 1201 . 

To summarize our analysis, the uP' , T^ behavior of 
the impurity electron selfenergy, which stems from the 
low-energy excitations and is at the heart of the Fermi 
liquid theory, is even quantitatively reproduced by the 
CTMA without spurious singularities. However, the lo- 
cation of this minimum at w = and the exact unitary 
value Adcr(O) of the impurity spectral density are not re- 
produced. The latter is crucial to avoid a non-causal 
behavior of the correlation part of the impurity selfen- 
ergy Tida{'^) and, hence, for an application as an impu- 
rity solver within DMFT^ . As will be discussed below in 
more detail, both of these failures can be attributed to an 
incorrect treatment of non-singular potential scattering 
processes at high energies, lo ^ Tk- We will propose a 
corresponding correction scheme in the next subsection. 

We conclude the present subsection by considering 
briefly another test case of Fermi liquid behavior, the 
empty orbital regime. In this case the d— electron den- 
sity of states consists of only one broadened, unoccupied 
single-particle resonance at a renormalized impurity level 
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FIG. 9: Logio-logio plot of Im Edo- — iO) — Im Y:da [i^o — iO) 
versus frequency, {uj — luo)/Tk, for T = 0.16 Tk- <jJo is 
the position of the minimum of Im T^da (i^ — *0) in Fig. |H| 
The dashed line is the fit to the low-frequency quadratic 
behavior of ImStjcr(w — iO) and represents the function 
y = 0.02440(a; - uof /T]^, see the text. 




T/Tk 



FIG. 10: Logio-logio plot of the minimum value of the 
imaginary part of the selfenergy, Im Y^da = <^min {T) , T) — 
lmT,da{uJ ^ ujo,T = O.OlTx) versus T/Tk- L>Jmin{T) is the 
position of the minimum for a given temperature T, see Fig. 
|S| Quadratic in T behavior is clearly visible for T < Tk- The 
dashed line is given hy y — 0.013 D (ttT/Tk) ■ 



Ed > far above the Fermi 
representation, the empty 
trivial strongly correlated 
ator constraint Q ~ \. Fif 
CTMA spectral density fo 
NCA is well known to fail 
a spurious, singular peak 
the X-ray-like divergences 



energy. In the pseudoparticle 
-impurity case is still a non- 
problem because of the oper- 
l- ini shows the NCA and the 
r ea/D = -hO.81, V/D = 0.2. 
badly in this case, producing 
at w = 0, which arises from 
of the pseudoparticle Green's 



functions. In the CTMA solution, the vertex corrections 
tend to cancel the infrared peak of the f — h bubble. Pre- 
sumably, the wiggles visible in the CTMA spectrum are 
due to numerical imprecision, but may also be due to 
a systematically imperfect cancellation of infrared diver- 
gent terms. In any case, the CTMA does not produce 
any definite peak structure at w = 0, significantly im- 
proving the description of the Fermi liquid behavior in 
the empty orbital regime. 




FIG. 11: CTMA (solid line) and NCA (dotted line) resuh 
for the local spectral density in the empty impurity regime, 
ed/D = -h0.81, V/D = 0.2 for T = 1.0 x 10"*D. The re- 
sulting occupation per spin is rido- = 0.037. The insets show 
the CTMA (left panel) and the NCA (right panel) spectral 
functions over the complete band width. 



B. Effective potential scattering correction 

Within the pseudoparticle technique, it is a non-trivial 
task to separate the single-particle hybridization part iV , 
of the total impurity selfenergy, from the interaction con- 
tribution, Y,da{^)^ since the hybridization is an interac- 
tion term in this representation, Eq. Hence, while 
the auxiliary particle method is designed for a systemati- 
cal treatment of the low-energy spin scattering processes, 
it is difficult to accurately calculate the non-singular po- 
tential scattering part of the total impurity selfenergy 



= iV 



ep - - GdA^r' (22) 
-Srf^M, (23) 



which involves hybridization processes at high energies 
of order e^. Technically speaking, the u; — value of 
/mE^'^*(ijj) is influenced by both, the real and the imag- 
inary part of G(jo-(0), and hence by the high-energy fea- 
tures of the spectrum through the Kramcrs-Kroenig re- 
lation. Obtaining the precise values of the real and imag- 
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inary parts of E^'^* (0 — iO) would therefore require calcu- 
lating the high-energy features of the local spectrum to a 
precision better than ~ Tk- Clearly, this is a formidable 
task, both with respect to numerical precision and to dia- 
grammatical systematics: Since any potential scattering 
term gives a non-singular, energy independent contribu- 
tion to S^'^'(O), it is unlikely that a class of principal 
diagrams can be identified that reproduces the correct 
value. On the other hand, the class of CTMA diagrams 
does describe the correct w^, behavior, reflecting the 
correct low-energy many-body dynamics. 
Based on these considerations, we propose a simple, phe- 
nomenological scheme to incorporate the correct poten- 
tial scattering contributions. It amounts to adding an 
appropriate frequency and temperature independent con- 
stant Ae to the real part of 'S*£^{uj). It has the effect 
of shifting the zero of the frequency scale in all quan- 
tities, and in particular in ImSc;cr(w), by virtue of the 
self-consistency42, Therefore, Ae can be chosen such that 
at T = the minimum of ImEdcri^-^ — *0) is obtained 
at cj = in accordance with Fermi liquid behavior, Eq. 
()19|l . In the CTMA solution the position of the minimum 
of ImSdai'-^) does not significantly change with temper- 
ature for T < 0.2Tk, see Fig. |S1 as expected from the 
Fermi liquid behavior, Eq. (|19|l . Hence, we have deter- 
mined Ae from the solution at T = O.IQTk to fulfill the 
Fermi liquid condition above. The results for the impu- 
rity spectral function corrected in this way is displayed 
for various T in Fig. El It shows accurate agreement 
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tures the Fermi liquid dynamics of the problem, miss- 
ing only part of the potential scattering contributions. 
The corresponding imaginary part of the total impurity 
selfenergy, Eq. (|22l) . is shown in Fig. E| Again, the 
Fermi liquid behavior, Eq. (|19|1 . is well obeyed. The min- 
imum value of Im'S*J'^{uj — iO) approaches for T ^ 
the value T^ff « 0.139 D instead of the exact limit 
r = 0.2 D. As discussed above, we attribute this to an 
inaccurate treatment of single-particle hybridization pro- 
cesses within CTMA. Note however, that for the DMFT 
algorithm^i only the interaction part 'Sdai^) of the self- 
energy is important. In the auxiliary particle method it 
is obtained from the impurity Green function by the sub- 



traction Sdcr(w — iO) = E 



lOJ- 



0) - iV 



where E'°* is 



given by Eq. H22|) , and its imaginary part remains strictly 
non-negative. 
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FIG. 13: Real and imaginary part of the total impurity 
selfenergy Y:^^J{lj — iO), corrected by the effective potential 
scattering method. The curvature of ImE^^*(Lj — zO) at the 
Fermi energy does not change considerably from T = O.IGTk 
down to T = O.OITk- The small bumb in the curve for 
T — O.OlT/f at positive uj is attributed to numerical inac- 
curacy. At Lj = 0,T ^ 0, ImE;,°'(0 - iO) assumes the ef- 
fective hybridization Fe// ~ 0.139D, somewhat smaller than 
the exact value F = 0.2D. Fe// is to be subtracted from 
ImE^'^*(a; — iO) in order to obtain the interaction part of the 
local selfenergy. 



FIG. 12: CTMA spectral functions with a constant potential 
scattering term added to the real part of T,scr(io — iO) according 
to the effective potential scattering method, see text. No low- 
energy singularity occurs, and the unitarity limit (arrow) is 
accurately fulfilled at the lowest temperatures. 

with the unitarity limit, even though this was not di- 
rectly implied by our adjustment procedure. This can be 
seen as a further indication that CTMA correctly cap- 



V. TWO-CHANNEL KONDO BEHAVIOR 

To complete the discussion of dynamical quantities we 
calculate the local spectral function of the two-channel 
{N — 2, M — 2) Anderson model. Here the low- 
temperature fixed point is of a distinct non-Fermi liq- 
uid natureiSi, involving a non-vanishing zero-point en- 
tropy, S{0) = ks ln-\/2, and a logarithmic divergence of 
the static susceptibility, x{T) oc — ln(T/Tx), signahng 
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a non-degenerate ground state and over-screening of the 
local spin, respectiveljsiJaSP . 

For the two-channel Kondo (2CK) model, the effective 
low-energy model of the two-channel SIAM, it has been 
shown using conformal field theory^ ^ that the local spec- 
trum has a cusp at the Fermi level, A^^^ {lj)~A^'~^^ (0) cx 
— Icl'I"'^/^. The weight of this power law becomes asymmet- 
rical for w > and w < 0, when the particle- hole symme- 
try is broken, e.g., by an additional potential scattering 
term. This weight asymmetry is analogous to the shift of 
the Kondo resonance id in the single-channel case. Ex- 
trapolating the Fermi liquid results of Appendix ^ to 
the two-channel SIAM, the weight asymmetry may be 
expected to be of OiTx/ed)- Very recently, the aux- 
iliary particle threshold exponents for the two-channel 
Anderson model have been calculated using the Bethe 
ansatziS It was shown that the exponents, like in the 
single-channel case, are functions of the local valence. 

The Bethe ansatz solution shows that the 2CK ground 
state involves intricate correlations between both con- 
duction channels and the local spin. Thus, one would 
expect that in a diagrammatic treatment three-particle 
correlation functions are needed, and that CTMA, which 
involves only two-particle T-matrices, is not able to cap- 
ture the correct 2CK ground state. Indeed, in the multi- 
channel (M > 2) case NCA as well as CTMA give in- 
correct, valence independent auxiliary particle threshold 
exponents as given in Section Surprisingly, how- 

ever, NCA correctly reproduces qualitatively the lead- 
ing low-energy singularities of physical quantities like the 
susceptibilitji^ or the local density of statesi^ In Fig. 1141 
we show the CTMA solutions for the impurity spectral 
function of the two-channel SIAM in comparison to the 
NCA result, both showing a cusp. The cusp of the 

NCA curve has a strong weight asymmetry, which is pre- 
sumably an overestimation, like in the Fermi liquid case 
(Fig. inj. This is significantly improved by the CTMA 
solutions. However, as mentioned above, we believe that 
a systematical description of the 2CK behavior would re- 
quire an extension of the CTMA to include three-particle 
correlation functions. 



VI. SUMMARY 

To conclude, we have extended the analysis of the Con- 
serving T-Matrix Approximation to thermodynamical 
and dynamical properties. It had been demonstrated ear- 
lier that for the single-channel Anderson impurity model 
the CTMA captures the correct spin-screened Fermi liq- 
uid ground state on the level of the auxiliary particle 
dynamics, signaled by the Fermi liquid exponents of the 
auxiliary particle propagators for all fillings^. In the 
present work we have shown that the CTMA also de- 
scribes the Fermi liquid strong coupling behavior of phys- 
ical quantities correctly. It is, thus, the first diagram- 
matic method that captures both the low-energy Fermi 
liquid behavior and the high-temperature properties of 
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FIG. 14: CTMA results for the impurity spectral density of 
the two-channel SIAM, with ed = -0.81 D , F = 0.2 D at var- 
ious temperatures. The NCA result is shown for comparison. 
The |<jj|^''^ cusp develops at the Fermi level, with a weight 
asymmetry due to potential scattering. The inset shows the 
spectrum over a wider frequency range. 



quantum impurity problems of the single-channel An- 
derson model type on the same footing^. In particular 
the CTMA describes the static spin susceptibility in the 
Kondo regime correctly for all temperatures. The Wilson 
number obtained within CTMA is in remarkably good 
agreement with the exact one. We also showed, that the 
physical d-electron spectral density in the single-channel 
case {N = 2, M = 1) is able to mend most of the de- 
ficiencies of the NCA in the Kondo and empty orbital 
regime of the model. Especially, no spurious infrared 
singularities occur. The deviation of the CTMA solution 
from the exact unitarity limiting value for the spectral 
density at the Fermi level could be traced back to an 
insufficiently accurate treatment of high energy poten- 
tial scattering processes. To correct this deficiency, we 
have proposed a phenomenological method by adding an 
appropriate, effective potential scattering term 5t to the 
real part of the impurity selfenergy, and defining an effec- 
tive single-particle hybridization rate ^eff- As a result, 
all essential Fermi liquid properties are fulfilled without 
spurious non-causal behavior. 

Finally, we comment on the prospects for future ap- 
plications of the CTMA. At the expense of being nu- 
merically involved, the CTMA combines two features 
which are non-trivial to fulfill by one single technique: 
flexibility and systematic treatment of the low-energy 
excitations without ad hoc assumptions about the na- 
ture of the ground state. These features may make the 
CTMA an attractive method for more complicated im- 
purity problems, such as (1) the selfconsistent quantum 
impurity problem that arises within the DMFT scheme^ ; 
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(2) quantum impurities with complex orbital structure; 
these arise also in cluster and cellular extensions of the 
DMFTiS^iSSiS; (3) quantum impurity problems which 
may exhibit a Fermi liquid instability. As a diagram- 
matic method, the CTMA is readily generalized for an 
arbitrary, energy dependent conduction electron DOS 
arising from the selfconsistent DMFT scheme. In ad- 
dition, the case of finite Coulomb repulsion U must be 
considered, in order to account for the upper Hubbard 
band and, e.g. to describe the metal-insulator transi- 
tion in the Hubbard model near half filling. It requires 
treating the bare charge fluctuation processes involving 
the empty and the doubly occupied impurity state in a 
symmetrical way. On an NCA-like level (Symmetrized- 
U NCA, SUNCA), this has been implemented in Ref. 
[ssf. see also [s^, and the corresponding Symmetrized, 
finite-U CTMA (SUCTMA) equations are reproduced in 
Ref. l28l|. The SUCTMA essentially amounts to cal- 
culating, in addition to CTMA, the ladder diagrams of 
heavy bosons representing the doubly occupied impurity 
state. Hence, the evaluations appear numerically man- 
ageable. Treating a more complex impurity orbital struc- 
ture requires introducing an individual auxiliary boson 
or fermion field for each charge and spin configuration 
of the impurity. Multi-orbital impurities have recently 
been treated using NCA-like approximationgSfliSi. The 
fact that the number of impurity configurations increases 
roughly exponentially with the number of orbitals will, 
however, limit the CTMA und SUCTMA to problems 
with not a too large number of local orbitals. On the 
other hand, because of the systematical treatment of low- 
energy excitations, the CTMA should at least be sensitive 
to instabilities of the Fermi liquid ground state due to, 
e.g., a quantum critical point. Future developments are 
planned to explore the possibility of these applications. 
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APPENDIX A: CTMA EQUATIONS 



In the following we present a compilation of the CTMA 
equations which determine the imaginary part of the im- 
purity propagator Gd, the auxiliary particle selfenergies 
T,f and Eh and the basic building block of the CTMA, 
i.e. the analytically continued 4-point vertices T*^"^^) and 
y(cfc)^ Due to the fact that the first diagram in FigQ] 
does not lead to a proper vertex contribution, we start 
the summation with the two-rung diagram. In order to 
include the proper channel and spin summations with 
enter even-rung and odd-rung diagrams differently, it is 
necessary to solve for the ladders with alternating signs as 
well. As already mentioned in the main text the two- rung 
diagrams, that is the inhomogeneous parts /(^^) and 
of Eqs. IAll and fA2l lead to selfenergy contributions which 
are already included in the NCA. Hence, in order to avoid 
over-counting these non-skeletons have to be subtracted 
before calculating the selfenergies. The resulting equa- 
tions are diagrammatically depicted in Fig^| We label 
the external frequencies of the 4-point vertices V^'^^^ and 
y'-'^''^ such, that the first (second) argument denotes the 
in- (out-)going frequency of the pseudo particle propa- 
gator and the third frequency labels the center of mass- 
frequency propagating through the ladder, see Fig. ^1 
At the Matsubara frequencies, the vertex functions are 
therefore given by 



V 



(Al) 



r 2 1 
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and 



^.1 



r 2 1 



After analytical continuation to the real axis the T- equations of the second kind: 
matrices obey the following linear Fredholm integral 



and 



f (If 



_ .(c6)(±) 



de 



r de 

^TtN{0) J TT 

r 



/(e - rt)GfAu; + e - 0)Gb^(£)AO,,(e - n) Gj^J + £ - 17). 



(A3) 



(A4) 



In order to simplify the expressions for the selfenergies with 
S/^CT and Sfc it proves useful to introduce 



and 



iV+ 1 
2 

iV+ 1 



( + ) ^ ^~ ^ y (c/) (-) _ TV , 

yi'^f) (+) - v^W) (-) _ /(c/) 



E<~-)(c.) 



Mr^ / -/(-£), <.(e)Gfc^(^-e) 



-,(c6) 



- 1 y(cfc) (+) _^ ^'^ + ^ y(c6) (-) _ J(c6) 



2 

M - 1 



y(cfc) (+) _ 1 y(c6) (-) _ jicb) 



Then we obtain for the analytically continued advanced 
{ioj ^ Lu — iO = uj) selfenergies (Fig. [T?)ll : 



(c.) = E;r'(a;) + S(;^)(c.) + E)7(c.) (A5) 



{u) = E-;-(c.) + s£^)(c.) + E--(c.) (A6) 



<cb), 



EfPiu;) = M i'^f{e-u)Al{e-u:) 



de 



_Afr / ^ f^fie-u:)fie'-c.) 

J TT J TT 

X yl°(u; - e)G6(e)7r7V(0)T2*"^^(£, e', £ + e' - w) 

X^O(c.-£')Gfc(£'), 
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where the vertex functions follow from T^^'^ and T^'^ 
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FIG. 15: Diagrammatic representation of the integral equa- 
tions comprising the T-matrices that enter in the expressions 
for the auxiliary particle selfenergies, see the text for details. 





FIG. 17: Diagrammatic representation of the equation for 
the local d-electron Green function. The first diagram is the 
NCA contribution. As described in the text, because of the 
analytical structure of the Adsigma it is necessary to multiply 
with a rung at each end of the vertex functions V^-^^^ , V^^'^^ , as 
shown in the last two diagrams on the right-hand side. Hence 
the sum of these diagrams contains at lease 4 rungs, i.e. the 
three-rung term must be added separately (second diagram 
on the right-hand side.) 



after performing the frequency summations and analyti- 
cal continuation of the external frequencies. In the con- 
tour integrals around conduction line cuts we have used 
AOH = (l7r)ImG2^,(o.-zO)/iV(0). 
Once the pseudo particle Green functions have been self- 
consistently calculated, physical properties can be deter- 
mined. In order to construct the local Green function 
from the basic building block of our theory, the 4-point 
vertices T'-'^f^ and T^^''', we note that the Bethe-Salpeter 
equations above yield either fully advanced or fully re- 
tarded T-matrices, 



T{iLU uj + iO, iLu' 



(A7) 

Lo' + iO,in n + iO) 



or 



FIG. 16: Auxiliary particle selfenergies E/ and E;,: The 
first diagram is the NCA contribution, where however fully 
CTMA-renormalized Green functions have to be used. The 
second and third diagram constitute the vertex corrections 
from T^f"^ and r''"'. 

and 

X ^°(e - uj)Gf (e) TTN{0)T^''^\e, e', e + e' - oj) 
x^O(£'-c.)G/(e'), 



T'^^'^{Lu,Lu',n)= (A8) 
T{iuj —* uj — iO, iLo' ^ Cij' — iO, ifi — > — iO) 



The local spectral density however requires the deter- 
mination of T^^^{ijo,uj' Therefore, its calculation is 
more involved. We construct T^^^(a;, w', fi) by adding 
two rungs to each of the diagrams in T^"^/^ and T'^'^''^ 
This then turns for example the two-rung diagram into 
a four-rung ladder. In order to include all contributing 
diagrams, it is then necessary to employ the T-matrices 
V^^^ and V^^^ instead of T^^f^ and T^^*) and, in ad- 
dition, the three-rung diagram has to be added, see Fig. 
1171 The equation for the advanced spectral function after 
projection onto the physical subspace then follows as 
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with 



AA(0) 
7r2r 



de- 



-/3e /-oo 



X Im{G/^(ei) Gf„{e2) Gb{ei + £2 - ?)} • Im {G/^(e - cj) Gb{ei - uj) Gb{e2 - t^)} 

/3c /-oo 



00 '^V(^) 



/ /(ei - i) /(e2 - - ^i) Al„^(i - e^) 

V_l.,(ei,e2,?)]} 



X Im{G/cr(e- w)G6^(ei -a;)Gh,,(e2 -w)} 

„-/3e /-oo 



00 '^^ f{-^) 



de, / de2 /(ei - ?) /(e2 - e)A'aa^{e, - i) A^^Je^ - e) 



X Im{Gfcp(ei)Gfcp(e2) 



,Af- 1 



F+^^,^,(ei,e2,e~) + 



Af + 1 



Kf:'.,,(ei,e2,e)]} 



X Im{Gfcp(e + tj)G/^(ei + w) G/^ (e2 + t^)} 



(A9) 



(AlO) 



APPENDIX B: FERMI LIQUID RELATIONS 

In this appendix we compile some exact Fermi hquid 
properties for the d-electron Green's function {N = 2, 
M = 1), 



Gda{u}±lQ) 



1 



(Bl) 



The interaction contribution Sd(T(w±iO) to the selfenergy 
obeys the Fermi hquid relations 

ImEd,(t^±zO) = ar^^^^^§^ , c^,T<Tk(B2) 

du;—^GdA^)=0 , T = (B3) 

These imply the Friedel sum rule, UdcriO) ~ S^, with 
ndcriO) the impurity occupation at T = and So- the 
scattering phase at the Fermi energy ep, 



cot S„ 



ReGrf^(O) 
ImGd^(0 + iO) 



(B4) 



The impurity spectral density at the Fermi energy w = 0, 
T = follows from Eqs. ||B1J, (jB^I, as 

A.. = ilmG..(0-0) = H^^^!(^. (B5) 

The Fermi liquid relations also determine the width 
r, position 6(1, and spectral weight z of the Kondo 



resonance. Using the quasiparticle weight z — {1 ~ 
dJ^dc7/d^\ui=o)~^, the Green's function can be written 
for \u;\, T<Tk, 



GdAi^±iO) = 



uj — id i iT 



(B6) 



with id = z[eF + ^d + R.el](jcr(0)], F = zT. z can be ex- 
pressed in terms of the strong coupling parameters of the 
Anderson model by equating the T = spin susceptibil- 
ity, Eq. (|12l) , with the Pauli susceptibility of a Fermi liq- 
uid of quasiparticles with weight z, x(0) = {gusY Ada / 



and hence, 



4r* sin'(7rnd^) 
ttTW 



r = — -- sin^(7rn<i^)T^ 
-kW 



(B7) 



(B8) 



Combining Eqs. HB4() . (|B6|I . HB8|I the position of the 
Kondo resonance relative to ep is given by 



= — F77 sin(27md^)T^ 



(B9) 



The exact prefactor a of the quadratic behavior of 
luYYjda{^) has been calculated for the Anderson model 
using direct perturbation theory to infinite order in 
the on-site repulsion [/^SiSiSi. One obtains at T = 
0, 52lmEd,(tj - zO)/ac^2 ^ ^|^(0,0)|2[A<i,(0)]3, where 
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7(0, 0) is the full, local 2-electron vertex function at the 
Fermi energy and at T = 0. Using a Ward identit3iS, 

it can be related to the quasiparticle density of states i^'^W'^ (i? — 1)^ \td\ 

Ad,(0)/z and expressed as 7(0, 0) = (i? - l)/[zAd,(0)], « = ^^3/2+20 sin^f^^,) ~d" ' 

where i? = 2 is the Wilson ratio. Hence, taking into ac- 
count Eq. we obtain for the prefactor in Eq. (|B2p . 
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